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ABSTRACT 

We reinvestigate the massive stellar content of the irregular dwarf galaxy NGC 3125 
(Tol 3) using VLT/FORS1 imaging and spectroscopy, plus archival VLT/ISAAC, 
HST/FOC and HST/STIS datasets. FORS1 narrow-band imaging confirms that the 
NGC 3125 A and B knots represent the primary sites of Wolf-Rayet (WR) stars, whilst 
HST imaging reveals that both regions host two clusters. Both clusters within region 
A host WR stars (Al and A2), for which the optically fainter cluster A2 is heav- 
ily reddened. It is not clear which cluster within region B hosts WR stars. Nebular 
properties are in good agreement with previous studies and infer a LMC-like mctal- 
licity of log(0/H)+12~8.3. LMC template mid-type WN and early-type WC spectra 
are matched to the observed blue and red WR bumps of Al and B, permitting the 
contribution of WC stars to the blue bump to be quantified. From our FORS1 spec- 
troscopy we obtain N(WN5-6 WC4)=105:20, ~55:0, 40:20 for clusters Al, A2 and 
Bl + 2, respectively. Our results are a factor of ^3 lower than previously reported 
by optical studies as a result of a lower Ha/H/3 derived interstellar reddening. Using 
Starburst99 theoretical energy distributions to estimate O star populations for each 
cluster, we find N(WR)/N(O)=0.2 for Al and 0.1 for A2 and the clusters within region 
B. From Ha narrow-band imaging, the O star content of the Giant Hn regions A and 
B is found to be a factor of 5-10 times higher than that derived spectroscopically for 
the UV/optically bright clusters, suggesting that NGC 3125 hosts optically obscured 
young massive clusters, further supported by VLT/ISAAC K band imaging. Archival 
HST/STIS UV spectroscopy confirms the low interstellar reddening towards Al, for 
which we have determined an SMC extinction law for NGC 3125, in preference to an 
LMC or starburst law. We obtain N(WN5-6)=110 from the slit loss corrected Hen 
A1640 line flux. This is in excellent agreement with optical results, although it is a 
factor of 35 times lower than that inferred from the same dataset by IChandar et alJ 
(2004j). The discrepancy is due to an anomalously high interstellar reddening derived 
from their use of the generic starburst extinction law. Highly discrepant stellar pop- 
ulations may result in spatially resolved star forming regions from UV and optical 
studies through use of different extinction laws. 
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1 INTRODUCTION 

In order to understand galaxy formation and evolution we 
need to be able to ac curately map the st ar formation his- 
tory of the universe jMadau et alj FL996l. Central to this 
topic are starburst galaxies, a class of object which display 
characteristics associated with massive, violent bursts of star 
formation. In the local Universe only a handful of starburst 

* Based on observations made with ESO telescopes at the 
Paranal Observatory under programme ID 074.B-0108 and with 
archival ESO VLT and NASA/ESA Hubble Space Telescope data, 
obtained from the ESO/ST-ECF Science Archive Facility, 
f E-mail: l.hadfield@shcf.ac.uk 



galaxies are res ponsible for a qu arter of the entire high-mass 
star formation (Heckman 1993). We need to be able to un- 
derstand nearby starbursts if we are to interpret the ob- 
servations of dis tant star forming r egions e.g. Lyman break 
galaxies fLBGs. Isteidel et alJll99d> . 

A subset of these galaxies are called "Wolf-Rayet galax- 
ies", since their integrated spectra display the broad emis- 
sion signatures associated with Wolf-Rayet (WR) stars; the 
highly evolved descendants of the most massive O stars. 
Spanning a wide variety of morphological types, WR galax- 
ies are observed in a wide variety of environments on a lo- 
cal (D<100Mpc) scale and out to high redshift. WR stars 
are exclusively associated with young stellar populations 
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(~5Myr), so WR galaxies represent an excellent diagnostic 
of recent star formation. 

Recently, IShaolev et alJ l)2003l) demonstrated that the 
composite spectrum of z~3 LBGs displayed broad Hell 
A1640 emission consistent with the presence of WR stars. 
The recent ultraviolet HST/STIS spectra l survey of lo- 
cal starburst galaxies (IChandar et al.ll2004l) revealed weak 
Hell A1640 emission in most cases. Of the 18 galaxies in- 
cluded in their surv ey, the super star clu ster NGC 3125-1 
(alias NGC 3125-A in lvacca fc Contil992]) showed the most 
prominent Hell A1640 emission. 

NGC 3125 (Tol 3) is a nearby (D=11.5Mpc 
ISchaerer et alJ I1999T) blue compact dwarf galaxy. Ob- 
servations have shown that the galaxy is dominated by a 
central starburst region which consists of two main emis- 
sion knots, NGC 312 5-A and -B. From UV spectroscopy, 
IChandar et alJ i2004l) estimate a WR population of ~ 5000 
and N(WR)/N(0)> 1 for NGC 3125-A. The latter is 
completely unexpected for the LMC-like metallicity of 
NGC 3125. In contrast, optical studies of NGC 3125-A infer 
a WR population of only ~ 500 and N(WR)/N(0)~ 0.1, an 
order of magnitude lower. WR populations in other nearby 
galaxies common to optical and UV surveys are found to 
be consistent to within a factor of two. If NGC 3125-A is a 
local analogue for LBGs it is necessary to reconcile optical 
and UV line techniques for this galaxy. 

Here we reinvestigate the WR population of NGC 3125 
using new VLT/FORS1 imaging and spectroscopy, supple- 
mented by archival VLT/ISAAC imaging and HST imaging 
and spectroscopy. This paper is organised as follows: VLT 
and HST observations of NGC 3125 are discussed in Sec- 
tion [5] Section [3] describes nebular properties derived for 
the two WR clusters using optical diagnostics. In Section 01 
WR and O star populations for each cluster are estimated 
using VLT/FORS1 imaging and spectroscopy. In addition, 
the WR population of NGC 3125-A is estimated from UV 
HST/STIS spectroscopy. Finally, we draw our conclusions 
in Section [S] 



2 OBSERVATIONS AND DATA REDUCTION 

We have observed NGC 3125 with the ESO Very Large Tele- 
scope UT2 (Kueyen) and Focal Reduced/Low dispersion 
Spectrograph #1 (FORS1), a 2048 x 2046 pixel Tektronix 
detector. Observations were made using the high resolution 
collimator which covers a 3.4' x 3.4' field of view, with a plate 
scale of 0.l"pixel _1 . Photometric observations of NGC 3125 
were acquired during November 2004 and January 2005, 
with spectroscopic data following in February 2005. Details 
of the observations, including the DIMM seeing can be found 
in Table 

To supplement VLT/FORS1 observations we have re- 
trieved VLT/ISAAC imaging and HST imaging and spec- 
troscopy of NGC 3125 from the ESO/ST-ECF archive. 

2.1 VLT Imaging 

FORS1 was used to obtain narrow-band images centred on 
A4684 and A4781 (FWHM=66A and 68A respectively). The 
A4684 filter coincides with the strong Nm (A4640A), Cm 
(A4650A), and Hen (A4686A) WR emission lines, whereas 



Table 1. VLT/FORS1 observation log for NGC 3125 



Date 


Observation Exposure 


DIMM seeing 






(sec) 


(") 


Imaging 


2004-Nov-08 


A4684 


60, 900 


0.42-0.43 




A4781 


60, 900 


0.42,0.38 


2004-Nov-17 


B 


10, 60, 600 


0.58-0.59 




A6665 


10, 60, 600 


0.45-0.57 


2005- Jan- 17 


A6563 


10, 60, 600 


0.67-0.79 


Spectroscopy 


2005-Jan-01 


300V 


4x 600 


0.60 



the A4781 samples a wavelength region relatively free from 
emission, providing a measure of the continuum level. In ad- 
dition, narrow-band on- and off-Ha images (A 6563, 6665A, 
FWHM=61,65A) were acquired along with broad-band B 
images. 

Images were reduced following standard reduction pro- 
cedures (i.e debiased, flat field corrected and cosmic ray 
cleaned) using iraf and STARLINK packages. We present con- 
tinuum, net Ha images and net A4684 images obtained with 
FORS1 in Figure 

The optical app earance of NGC 312 5 is that of an 
amorphous elliptical (ISchaerer et al.lll999l) . dominated by 
a bright central starburst. In Figure 0^a) we present a B 
band image of the central starburst region taken in excel- 
lent seeing conditions (~0.6"). The continu um emission is 
domi nated by two emission knots A and B iVacca fc Contil 
1992). A corresponds to the slightly brighter knot and is lo- 
cated ~10" to the NW of B. In the image, knot A appears 
to be partially resolved into two components whereas knot 
B does not. Several fainter emission knots appear to link the 
two dominant regions. 

Figure □» shows the VLT/FORS1 net Ha image of 
the central starburst. Nebular emission is concentrated on 
the two principal knots with region A extending 300pc and 
region B 200pc. This large size and Ha luminosities of 
3.2xl0 40 and 2.0xl0 40 ergs" 1 for A and B, respectively (see 
Section 14.2.21 . indicate that A and B are typical of extra- 
gala ctic Giant Hu re gions which host multiple stellar clus- 
ters jKennicuttJll984h . 

We have searched for characteristic WR signatures 
NGC 3125 by subtracting the A4871 image from the A4684 
image, which is shown in Figure 0^c). Regions A and B are 
clearly identified as the primary hosts of WR stars. Region 
A contains a fainter emission component ~0.4"to the W of 
the main source. This fainter component contributes ~15% 
of the total A4684 emission. 

Photometry of knots A and B was performed using 
the aperture photometry package phot within iraf. In- 
strumental zero points were derived by observing LTT3864 
(B = 12.7) and GD108 (B = 13.3) spectrophotometric stan- 
dard stars. 

Magnitudes for source A include the fainter compo- 
nent to the W of the knot. We find niA4686 — 17.0±0.1 and 
m>478i = 17.2±0.1. Photometric uncertainties have been esti- 
mated from background variations, contamination by nearby 



How extreme are the Wolf-Rayet Clusters in NGC 3125? 3 



F0RS1/Bessell B F0RS1/net Halpha 




10:06:34.0 33.6 33.2 10:06:34.0 33.6 33.2 

Right ascension Right ascension 

Figure 1. 15" X 15" VLT/FORS1 and VLT/ISAAC archival images of NGC 3125. For our assumed distance 11.5Mpc, the physical 
region illustrated equates to 750 X 750 pc. a) Bcsscll B image showing the morphology of the central starburst. Regions A and B have 
been marked along with the slit position of spectroscopic observations, b) High contrast, net Ho image, c) Difference between A4684 and 
A4871 filters, showing HenA4686 / C in A4650 emission, d) Archival VLT/ISAAC K s image of NGC 3125. North is up and East is to the 
left on all images. 



emission sources and the aperture size. An extremely close 
emission knot to the NW of source B made aperture pho- 
tometry more difficult, reflected by larger photometric un- 
certainties, for which we derive mA4686 an d ni\478l magni- 
tudes of 17.5±0.2 and 17.7±0.2 respectively. 

To complement our imaging dataset we have re- 
trieved archival VLT/ISAAC images obtained on 18 April 
2002 under excellent seeing conditions of FWHM~0.4". 
These observations comprise broad-band K s imaging (5 
on source plus 2 off-source exposures of 5 s) taken us- 
ing the SW mode (Rockwell-Hawaii detector, plate-scale = 
0.148 arcsec pixel -1 ). Images were calibrated relative to the 
2MASS source 10063237-2956190 (K s =15.0mag) which was 
present in this 2.5'x2.5'field. 

In contrast to optical images, Figure^d) shows that in 
the near-IR, knot A is much brighter than B, with Ks=14.5 
and 16.1 mag, respectively. In addition to the optically dom- 
inant sources, there are two bright knots to the NW of 
B, with a combined magnitude of Ks=15.0mag. These are 



again more prominent than in optical images, suggesting the 
presence of partially obscured star forming regions. We will 
return to this in Section f4.2.2l 



2.2 VLT Spectroscopy 

Spectroscopic observations of NGC 3125 were undertaken on 
1 February 2005 using FORS1 and the high resolution colli- 
mator in seeing conditions of ~0.6". Spectra were acquired 
using a 0.8" slit and 300V grism centred along the two main 
emission knots (PA =-124°; see Figure |Jc)). The extracted 
spectra covered a wavelength range of 3300-8600A with a 
dispersion of 2. 6A pixel -1 and resolution of ~ 15A (as mea- 
sured from comparison arc lines). 

Data were prepared and processed using standard tech- 
niques and iraf and STARLINK packages i.e. bias subtracted, 
flat field corrected, extracted and flux / wavelength cali- 
brated. Care was taken during the extraction process to en- 
sure neighbouring emission knots did not contaminate the 
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Figure 2. VLT/FORS1 optical spectra of NGC3125-A and 
B. Spectra have been ve locity corrected (V r = 865km s —1 
lLauberts fc Valentiinll989D and WR emission features have been 
marked. 

aperture or the background subtraction. The spectrophoto- 
metry standard Feige 66 was observed in order to relatively 
flux calibrate the spectra. 

Absolute flux calibration was achieved by comparing 
synthetic-filter photometry (A c = 4684, FWHM = 66A) 
to A4684 photometry. We derive a slit correction factor of 
1.3±0.2 for knot A, which dominates the accuracy of our 
final flux calibrated spectra. An identical correction factor 
is found for knot B. We expect our absolute flux calibration 
to be correct to <~15%. Extracted spectra of NGC3125-A 
and B are shown in Figure [5] in which the principal WR 
emission features are identified, including weak N IV A4058. 

2.3 HST Imaging 

To compliment our FORS1 dataset we have retrieved 
archival HST images, including a STIS/LONG_PASS acqui- 
sition image (t = 40 s, A c = 7200A) from programme GO 
9036 (P.I. C. Leitherer) of region A and a FOC/F220W 1 im- 
age (t =500 s, A c = 2280A) from programme GO 4800 (P.I. 
P. Conti) of the central starburst region. 

The STIS/LONG_PASS acquisition image resolves re- 
gion A into two clusters, separated by ~ 0.5"or 25pc at 
the distance of NGC3125 (upper panel of Figure [3J. We 
designate these two components Al and A2, with Al corre- 
sponding to the brighter of the pair to the east. The com- 
ponents have FWHM of 0.17" and 0.13", and a flux ratio of 
Fai/Fa2 ~1.9 at 7200A, although an absolute flux calibra- 
tion was not possible. 

In the HST/FOC F220W image Al is bright whereas 
A2 is barely detected, with F A i(F220W)/Fa 2 (F220W) ~10 
(lower panel of Figure |3J. Assuming that Al and A2 pos- 
sess comparable intrinsic energy distributions - which is 
reasonable since both clusters host WR stars (Section 12. in 
- A2 must suffer from significantly higher extinction. 
Aperture photometry reveals F A i(F220W) = (l.l±0.2) x 
10 -15 ergs -1 cm -1 and F A2 (F220W) < 3x 10 -16 ergs -1 cm -1 . 

1 FOC/F220W image is pre-COSTAR. 
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Figure 3. l"x l" (50x50 pc) STIS/Long_Pass acquisition (top) 
and pre-COSTAR HST/FOC F220W (bottom) images of regions 
A. The two clusters Al and A2 are marked. North is up and East 
is to left. 



The HST/FOC image also resolves region B into two 
components (Figure|3J - Bl and B2, with the brighter com- 
ponent, Bl to the SE with Fbi(F220W)/F B 2(F220W)~1.2. 
It is not clear which cluster within region B hosts the 
WR population since these are spatially unresolved in our 
VLT/FORS1 imaging. High spatial resolution optical imag- 
ing is not yet available. It was not possible to perform aper- 
ture photometry on Bl and B2 individually in these pre- 
COSTAR images, we therefore derive a total UV flux of 
(1.6 ±0.2) x 10" 16 ergs" 1 cm" 1 for region B. 

2.4 HST Spectroscopy 

Archival HST/STIS spectroscopy of region NGC3125-A1 
llChandar et al jr2004T) was obtained in programme GO 9036 
(P.I. C. Leitherer) and comprised two exposures with the 
G140L grating (FUV-MAMA detector), of duration 1050 
and 2925 s, plus one exposure with the 230L grating (NUV- 
MAMA detector) using the 52x0.2"slit. The STIS MAMA 
detectors have a plate scale of 0.024" pixel -1 . Spectra were 
extracted over a 13 pixel (0.3") aperture and combined 
to provide a complete wavelength coverage of 1175-3100A, 
with a spectral resolution of ~3 pixels (1.8 A). 

The UV s pectr um of Al has been presented by 
IChandar et alJ l|2004l ) who noted the strong Hell A1640 
emission feature, for which we measure an equivalent width 
of 6.9±0.8A and FWHM~4.8±0.7A. Absolute flux calibra- 
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F0C/F220W 




Table 3. Summary of nebular properties of NGC 3125 clusters A 
and B. 



33.93 10:06:33.90 
Right ascension 

Figure 4. l"x 1" (50 x50pc) pre-COSTAR HST/FOC F220W 
(A c = 2280A) image of region B. North is up and East is to left. 



Table 2. Observed (Fa) and dereddened (Ia) nebula line fluxes of 
WR regions A and B within NGC 3125. Line ratios are normalised 
to H/3=100. We present observed and dereddened H/3 fluxes in the 
final row (erg s — 'cm" 2 ). Reddening corrections of Es_y-=0.24 
and 0.21 for A and B respectively, include a Galactic foreground 
reddening of 0.08 mag. 



•Vest 






A 




B 












(A) 


ID 


Fa 


Ia 


Fa 


Ia 


3727 


[On] 


83.5 


99.0 


109.2 


123.9 


4340 


H 7 


40.0 


44.7 


40.3 


43.3 


4363 


[Oin] 


3.9 


4.6 


2.8 


2.9 


4861 


H/3 


100 


100 


100 


100 


4959 


[O in] 


205.7 


203.2 


169.8 


162.3 


5007 


[Oni] 


649.3 


615.2 


524.2 


496.0 


6563 


Ha 


371.4 


282.9 


351.2 


279.2 


6583 


[N II] 


8.7 


6.6 


11.5 


9.2 


7330 


[On] 


3.7 


2.7 


3.6 


2.7 


4861 


H/3 


6.95x10" 


" 14 1.56x10' 


13 3.83x10" 


14 7.55X10" 14 



tion has been achieved by comparing our HST/STIS spectra 
with the F220W flux of Al. We derive spectral slit losses of 
-40±30%. 



3 NEBULAR ANALYSIS 

In the following section we will derive the nebular properties 
of NGC3125-A and B. Observed (Fa) and dereddened (Ia) 
nebular line fluxes with respect to H/3=100 are presented in 
Table H 

The nebular analysis was performed using the star- 
link package DIPSO with line fluxes being determined using 
the ELF (emission line fitting) routine. The emission line 
profiles were non-gaussian in the FORS1 spectra (compari- 
son arc lines displayed the same profile, apparently due to 
the use of the high-resolution collimator for spectroscopy) 
and were modelled using template profiles created from the 
data. [O ill] and [O n] emission lines were modelled using 
[O in] A5007 as a template profile, Ha, /3, 7 and [N 11] were 
modelled using the H/3 emission line as a template. 





A 


B 


Te(Oll) (K) 


12100±1000 


10 800±500 


Te(OlIl) (K) 


10 500±500 


9 800±500 


n e (cm" 3 ) 


140+ 


140+ 


O+/H 


(1.8±0.2) x 10~ 5 


(3.5±0.2) x 10~ 5 


2 +/H 


(1.9±0.1) x 10" 4 


(1.9±0.1) x 10~ 4 


log(0/H) + 12 


8.32±0.03 


8.35±0.03 



Electron densities were derived for [S 11] ratios of VC92. 



3.1 Interstellar reddening 

Estimates of the interstellar reddening have been made using 
the Balmer line ratios Ha: H/3: H7. Nearby [N 11] emission has 
been accounted for when measuring observed Ha fluxes. As- 
suming Case B recombination theory fo r electron densities 
of 10 2 cm -3 and a temperature of 10 4 K jHummer fc Storevl 
Il987l) we deduce average total Eb-v values of 0.24 and 0.21 
for A and B, respectively. 

Foreground Galactic red dening (Eb-v) to wards 
NGC 3125 of £ s _v-=0.08mag dSchleeel et alJ ll99gFl was 
accounted fo r using a standard Galactic extinction law 
l|SeatorJll979r) . The point-like appeara nce of the two regions 
on ground-based images suggest that a lCalzetti et al.l Jl994T) 
starburst obscur ation law is inappro priate. We therefore 
choose to use the lBouchet et all il985T) SMC extinction law 
to deredden our spectra (see Section f5. 211 . 

Underlying stellar Ha and H/3 absorption from early- 
type stars is estimated to be Wa ~2A. In our spectra, Ha 
and H/3 equivalent widths are found to be ~550A and ~100A 
for both regions. Propagating this correction through calcu- 
lations leads to an uncertainty in Eb-v of ±0.01 mag. 

From a comparison between the STIS/Long_Pass ac- 
quisition and FOC/F220W images (Figure EJ, cluster A2 
appears to suffer a much higher extinction than Al, with 
Fai/F A2 ~1.9 at 7200A and ~10 at 2280A. If we assume 
that Al and Al possess identical intrinsic flux distributions, 
we estimate that E™±\ / (A2)~ 0.5mag. 

iKunth fc SargentJ jl98ll) derived a total E B _v=0.40 
from H/3:H5:H7 intensity ratios 2 . Higher (and weaker) mem- 
bers of the Balmer series can be significantly affected by un- 
derlying stellar absorption, which was not accounted for in 
thei r analysis. 

IVacca fc Contl lll99l hereafter VC92) also studied the 
nebular properties of NGC 3125-A and B. Using Ha:H/3 line 
ratios they derived internal reddenings of E™^y=0.40 and 
0.64, respectively, significantly higher than those obtained 
here. A higher extinction for region B than region A is in- 
consistent with our UV and optical photometry, if we as- 
sume identical flux distributions for these regions, since B is 
brighter in the UV whilst A is brighter optically. 

VC92 based their determination on composite blue (A < 
5400A, photometric) and red (A >4500A, non-photometric) 



2 IKunth fc SargenJ 11981 
EB_v=0.28mag in their Table 1, the 
not given. It is not clear which value has been used in their 
subsequent analysis. 



quote 
Table 1, 



colour excess of 
origin of which is 
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spectra obtained 12 months apart. The red spectrum, which 
was scaled to the blue continuum, resulted in a H/3 flux 20% 
larger than that measured in the blue spectrum. Based upon 
their red dataset alone, an extinction of NGC 3125-A fully 
consistent w ith the present result would have been obtained 
iVaccall2005l private communication). 

3.2 Electron temperature, density & oxygen 
abundance 

Electron temperatures, T e , for regions A and B have been 
derived from the temperature diagnostics [O n] 3727/7325 
and [Om] (4959+5007) /4363 A for the line ratios presented 
in Table |21 Electron temperatures were calculated using the 
five-level atom calculator temden within iraf for a constant 
electron density of 100 cm -3 . 

Errors on T e [O n] were based on the 5% measurement 
error estimated for the [On] 7325Aline, while T e [Olll] un- 
certainties were based on the 10% formal error given for 
[O in] 4363A. 

Electron temperatures derived here agree to within the 
errors of those reported by VC92 who found 10 200 K for 
region A and 10 300 K for region B. 

Due to the low spectral resolution of our data 
[S n] 6716/6731A nebula lines are not resolved, preventing a 
direct estimate of the electron density n e . For completeness, 
we have estimated n e by combining our electron tempera - 
tures with [Sn] fluxes published bv IVacca fc ContH il992t) . 
From these we estimate an electron density of 140 cm -3 for 
both regions. Errors are not quoted for electron densities 
since published [S n] ratios did not give associated uncer- 
tainties. 

Since the oxygen content is used as a proxy of a galaxy's 
metallicity, we have derived the oxygen abundance for each 
region using [On] 3727A and [Om] 5007A nebular emis- 
sion lines and their associated electron temperatures (see 
Table 0J • 

Both regions, with log(0/H) + 12 = 8.32-8.35 have 
an oxygen content com parable to log(0/H) + 12 = 8.37 
jRussell fc Dopitai ri990) observed in the LMC. Our oxygen 
abund ances are in excell ent agreement with previous stud- 
ies of iKunth fc ' Sargentl and IVacca fc Contil who obtained 
log(0/H)+12 = 8.3-8.4. 

4 MASSIVE STAR POPULATION 

In this section we derive the massive star content of 
NGC 3125-A and B using optical and UV techniques. 

4.1 Estimating the number of WR stars 

4-1.1 Optical 

For regions A and B we detect br oad blue and red WR 
features, as previously reported bv ISchaerer et ail lll999l) . 
Gaussian line profiles have been fit to the Nm A4640/Cm 
A4650 blend, Hen A4686 and Civ A5808. Slit loss corrected 
emission line fluxes are presented in Table [I] 

In region A we detect strong He II A4686 emission plus 
Nm A4640 - Cm A4650, indicative of a predominantly late 
WN population. We also detect weak N IV A4058, as reported 



Table 4. Observed WR line properties and derived popula- 
tions of clusters in NGC 3125. Line fluxes (F\) are expressed 
in ergs - cm , deriv ed luminosities adop ted a distance to 
NGC 3125 of 11.5Mpc ISchaerer etal]|l99a) and are expressed 
in ergs -1 . 



Region 


Al 


A2 


Bl + 2 




0.24 


-0.58 


0.21 


F4686 


7.4xl0 -15 


~l.lxl0 -15 


4.0xl0 -15 


I4686 


1.7xl0 -14 


~6.8xl0 -15 


8.4xl0 -15 


L4868 


2.7xl0 3s 


~1.0xl0 38 


1.3xl0 38 


N(WN5-6)t 


150 


~55 


75 


F58O8 


2.0xl0 -15 




2.8xl0 -15 


I58O8 


3.8xl0 -15 




4.8xl0 -15 


L58O8 


6.0xl0 37 




7.6xl0 37 


N(WC4)t 


20 




25 


Fl640 


2.0xl0 -15 






Il640 


1.3xl0 -13 






Ll640 


2.1xl0 39 






N(WN5-6)t 


115 







' Numbers derived solely from o bserved line luminosities bas ed 
upon average line luminosities of lCrowther fc Hadfieldl J2006t) 



by IKunth fc Sargentl dl98ll). This suggests a dominant 
mid WN subtype (ISmith et alJll99rj ; ICrowther fc Hadfieldl 
l200rJ) . 

The presence of C IV A5808 emission suggests that 
a significant WC population is present in region A. WC 
star s are classified by the Cm A5696/Civ A5808 line ra- 
tio JSmith et alJll990all . Since C ill A5696 is very weak or 
absent, the dominant WC population is WC4-5. 

Similar WR features are observed for region B, although 
Hell A 4686 is somewhat weaker than that in region A. As 
such, we assume a dominant mid WN population whilst the 
red feature is again consistent with a early-type WC popu- 
lation. 

As a first estimate of the WR population in each re- 
gion we have derived numbers based solely on the observed 
Hen A4686 and C iv A5808 line lumi nosities, following the 
approach of ISchaerer fc Vaccal (11998^1 . We have attempted 
to derive WR populations for Al and A2 separately, since 
UV spectroscopy relates only to Al and A2 suffers a much 
higher extinction. From Figure^c) we estimate that the ob- 
served Hen A4686 flux ratio is 0.15:1 for A2:A1; whilst we 
assume only Al is responsible for the C IV A5808 emission. 
This assumption is based on the higher extinction of A2 
which suggests that it may be younger than Al and would 
not necessarily host a mixed WR population. 

If we assume that only WN5-6 stars contribute to the 
Hell A4686 line flux i.e. neglecting the WC contribution, 
we estimate that N^(WN5-6) = 150 based upon the (1.8 ± 
1.7) x 10 36 ergs -1 averag e Hell line luminosity of 15 LMC 
WN5-6 stars studied bv ICrowther fc Hadfieldl fc006h . For 
A2, assuming an internal reddening of E™^ v ~ 0.5 mag, we 
estimate that N A2 (WN5-6)~55. Similarly, Nai(WC4)=20 
based upon the average (3.3 ± 1.7) x 1 36 ergs - ^ C IV A5808 
line l uminosity of 7 LMC WC4 stars JCrowther fc Hadfieldl 
l2006h . 

Applying the same methodology to region B, we derive 
a WR population of 75 WN and 25 WC stars. 
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Figure 5. Dereddened (SMC extinction law), continuum subtracted, spectral comparison between the observed (solid) and generic 
(dashed-dotted) WR emission features for clusters Al and B. Observed spectra have been velocity corrected and corrected for 25% slit 
losses. Generic WC4 (dashed) and WN5— 6 (dotted) features are marked. 



Although WN stars will typically be the primary con- 
tributor to the Hell A4686 feature, there will be a con- 
tribution from the WC population if present. Therefore, 
to improve upon our initial estimate of the WN content 
we have estimated the WC contribution to the blue fea- 
ture by fitting generic LMC WN5-6 and WC4 spectra from 
ICrowther fc Hadfieldl i2006l) to the observed blue and red 
WR bumps. Each region shall now be discussed in turn. 

In Figure OJb) we compare the C IV A5808 profile of 
region Al with that expected from 20 LMC-like WC4 stars 
at a distance of NGC 3125. The good agreement confirms our 
initial estimate. We have then accounted for the contribution 
of 20 WC4 stars to the blue WR feature, adjusted for the 
contribution of A2 to the observed flux. 

In Figure 1^1 a) we compare the dereddened, continuum 
subtracted blue WR feature for cluster Al with that ex- 
pected for a mixed WR population of 20 WC4 and 105 
WN5-6 stars. The composite spectrum of our generic WR 
populations reproduces the observed morphology exception- 
ally well, except for spectral regions where nebular lines are 
expected (e.g. [Fein] 4658, [Anv] 4711). 

For region B, it was necessary to reduce the number of 
WC4 stars from 25 to 20 in order to match the observed 
Civ A5808 emission profile (Figure [jjd)). This revised WC 
population was used to construct the blue WR bump for 
region B. Figure Efc) shows that the WC contribution is 



highly significant, with 50% of the Hell A4686 and 90% of 
the Niii/Ciii A4640/50 line flux originating from WC stars. 
The WN5-6 population is therefore reduced to 40. WR pop- 
ulations derived using this method are included in Table |K| 

4.1.2 UV 

The WR population of cluster Al has been independently 
estimated from the dereddened, slit loss corrected Hell 
A1640 line luminosity. Adopting E_b-v=0.24 mag, as de- 
rived from Ha:H/3, we derive a dereddened line flux of 
F A1 (1640)=9.4xl0" 14 ergss" 1 cm" 2 based upon an SMC ex- 
tinction law. For completeness we have considered various 
extinction laws, and conclude in Section [5. 2l that an SMC- 
cxtinction law provides the closest match to the complete 
STIS UV spectral energy distibution. If we neglect the WC 
contribution to the A1640 line we estimate a stellar content 
of 115 WN5-6 stars, based upon the a verage Hen A1640 line 
lumin osity of 1.8xl0 37 ergs _1 from ICrowther fc Hadfieldl 
(2006). Application of an LMC extinction curve, which is 
less successful at reproducing the UV flux distribution of 
Al, requires a higher extinction of Es-v = 0.33 mag and a 
WN population of ~200 stars. 

Following the same approach as for the optical WR fea- 
tures, we have constructed a He II A1640 profile using generic 
LMC WN5-6 and WC4 spectra. Generic spectra are taken 
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Table 5. WR populations for clusters within NGC3125 A and B 
derived from fitting LMC template WR spectra. 



Region 


Diagnostic 


Al 


A2 


Bl + 2 


N(WN5-6) 


A4686 


105 


~55 


40 


N(WN5-6) 


A1640 


110 






N(WC4) 


A5808 


20 




20 


N(WR) 


A4686/A5808 


125 


~55 


60 



from ICrowther fc Hadfieldl l|200ffl and are based on low res- 
olution IUE/SWP data for 10 WN5-6 stars and medium 
resolution HST/FOS data for 6 WC4 stars. 

In Figure|^]we compare the SMC law dereddened, far- 
UV slit loss corrected HST/STIS spectrum of Al, plus 
the spectrum around Hell A1640 degraded to the 6A res- 
oluti on of IUE/LORES - with that for 110 WN5-6 stars 
from lCrowther fc Hadfieldl i2006t) adjusted to the distance of 
NGC 3125. Also illustrated is a synthetic Starburst99 spec- 
trum based upon LMC/SMC template stars for the 1250- 
1600A region, this will be further discussed in Section 14.21 
The WC contribution to the Hell A1640 flux is very minor, 
reducing the number of typical WN5-6 stars from 115 to 
110. This is in excellent agreement with the optically de- 
rived WN5-6 population (see Table |5Jf 

4.2 Estimating the number of O stars 

4-2.1 The O star population of individual clusters 

The O star content of Al has been directly estimated by 
comparing the slit loss corrected, dereddened HST/STIS 
spectrum of Al to the best fitting Starburst99 model (SB99, 
iLeitherer et al.1 ll99Sh . We assume a Kroupa IMF (0.1- 
IOOMq) with a turnover at 0.5Mq and an exponent of 2.3 
for the high mass interval. The Geneva high mass loss rate 
LMC metallicity model has b een adopted. For an e mpirical 
LMC/SMC template spectra (ILeitherer et al]l200l|) we esti- 
mate a burst age of 4Myr due to the prominent Siiv A1 400 
feature, see Figure |S| In contrast. IChandar et alJ <l2004f) es- 
timated 3±1 Myr. A cluster mass of 2. Ox 10 s Mq is required 
to match the UV continuum, indicating an O star content 
of ~550. The N(WR)/N(0) ratio for cluster Al is therefore 
~0.2, explaining the large equivalent width of Hen A1640. 

If clusters Al and A2 are coeval and using the flux ratios 
described in Section 12.11 we find that A2 is slightly more 
massive, with M=2.2x10 5 Mq. The O star population of A2 
is estimated at ~600, giving N(WR)/N(O)~0.1. 

For region B, the SB99 model was scaled to match 
the FOC(F220W) and VLT/FORS1 flux levels which in- 
dicates a combined mass of 1.6xlO 5 M . Given their rel- 
ative UV fluxes, as described in Section 12.31 and assum- 
ing each cluster has an identical E™^ v , we estimate that 
M(Bl)=9xl0 4 M Q and M(B2)=7x 1O 4 M . SB99 models 
predict that such bursts should host 250 and 200 O stars, 
respectively. The O star content derived for these clusters 
is 450 leading to N(WR)/N(O)=0.1, as in A2. A summary 
of the cluster O star content and N(WR)/N(0) ratios are 
presented in Tabic 

In Figure [7] we compare slit loss corrected VLT/FORS1 
and HST/STIS spectra with the combined, reddened, SB99 
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Figure 6. Comparison between the slit loss corrected, dered- 
dened (E£OT, =0 ,.08(GAL) + 0.16(SMC)) HST/STIS spectrum 
of Al (faint solid) with a 2.0x10 5 Mq SB99 synthetic spec- 
trum obtained using LMC/SMC template OB stars for 1250- 
1600A(solid), in which an age of 4Myr is indicated from the promi- 
nent Si IV wind profile. Also shown is the STIS spectrum around 
He IIA1640 degraded to 6A resolution (dotted) together with the 
A1640 e mission predicted from 110 generic LMC WN5— 6 stars 
(dashed, ICrowther fc Hadfieldlfeood) . 

models for clusters within regions A and B. Individual 
contributions for each cluster are indicated, together with 
transmission curves for the FOC/220W and STIS/LP fil- 
ters. In order for the SB99 flux distributions of each re- 
gion to match their counterpart VLT/FORS1 spectra, we 
estimate that E I ^ 1/ (A1)=0.15 and EgT v (B)=0.11, in ex- 
cellent agreement with those derived from nebular Ha/H/3 
and UV spectroscopy /photometry. Best agreement between 
STIS and FORS1 spectroscopy for region Al was achieved 
with a slit correction factor of 1.6 versus 1.4±0.3 obtain from 
HST/FOC imaging. 

4-2.2 The O star population of the Giant Hll regions A & 
B 

In addition to determining the O star content of the indi- 
vidual clusters from their UV/optical continuua, we have 
derived the number of O stars present in each giant H 11 re- 
gion using the net Ha flux measured from our VLT/FORS1 
narrow-band on- and off-Ha images. 

The observed net Ha fluxes of regions A and B are 
1.2 and 7.5 x 10~ 13 ergs" 1 cm -2 , respectively. For internal 
reddenings of 0.16 and 0.13 mag, these equate to Ha lu- 
minosities of 3.2 and 2.0xl0 40 ergs _1 . Measurements were 
made using apertures 7.0"and 6.0"in diameter for A and B, 
respectively. These correspond to a physical scale of 350 and 
300 pc. From slit spectroscopy, Ha luminosities are measured 
to be a a factor of ~3 lower. 

Accounting for the WR contribution to the ionising con- 
tinuum, the number of equivalent 07V stars, Noiv, can be 
expressed as 

»r Qo bS ~ NwnQ^) N + NwcQ^O C 
= QgTV 

where Qo bs is the observed Lyman continuum flux and 
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Figure 7. Left panel: Comparison between the observed VLT/FORS1 (A1+A2) and slit loss corrected HST/STIS spectroscopy of 
Al. Reddened flux distributions for Al (M=2.0x10 5 Mq - dashed) and A2 (M=2.2xlO 5 M0 - dotted) are shown together with the 
combined, A1+A2, flux distribution (solid). We adopt a foreground extinction of Es-v=0.08 mag, with internal (SMC law) extinctions 
for individual clusters indicated. The F220W photometry for Al and A2 and the relative STIS/Long_Pass flux ratio Fai/Fa2 ~1.9 plus 
the FOC/STIS filter transmission curves are also shown. Right panel: A similar comparison for clusters Bl and B2 assuming identical 
internal extinctions, plus the combined Bl+2 FOC/F220W flux. 



Qo N , Q^o , Qo ?V a re the average Lyma n continuum flux 
for each stellar type llVacca fe Contilll992l) . 

The O star population present in our regions will of 
course not be restricted to the 07V spectral type, but dis- 
tributed amongst the entire spectral class. We must there- 
fore account for the age of the population and the IMF when 
determining the total number of O stars (No)- As shown by 
IVaccal dl994f) . No is related to the number of equivalent 07V 
stars, No7v, by 



No = 



Noi 



where r/(t) is the IMF averaged ionising Lyman continuum 
luminosity for a stellar population, of an given age nor- 
malised to one equivalent 07V star. 

In recent years a numb er of papers (i.e. iMartins et al] 
2002:; ICrowther etUI [2002) have re-calibrated the spec- 
tral type-temperature relation for Galactic O stars, ac- 
counting for non-LTE and line-blanketing effects. Stud- 
ies have shown that a ~10% downwards revision in the 
effective temperature scale of O stars is required, such 
that log Qp for a typical G alactic 07V star decreases by 
0.2 dex jMartins et alJ 120 05). Additionally, recent studies 
of Magallenic Cloud O stars indicate 2-4k K higher tem- 
peratures than their G al actic counterparts {Massev et alJ 
120051: iHeap et alJ l2006t iMokiem etaO 120061) . Therefore, 
for this present application we adopt log Q07 y=48.9 for 
LMC metallicity stars (and would recommend 49.0 for 
SMC metallicities). For the WR contribution, we adopt 



log QcT C 



: 49.40 and loi 



49.75, based on the aver- 



se. 

age of 6 WC4 LMC stars Jcrowther et alj2002t and 9 WN5- 
6 LM C stars JCrowther fc SmitrJl996tlCrowther fc Dessard 



The parameter r\ has been evaluated using t he in - 
stantaneous starburst models of ISchaerer fc Vaccal <ll99St) 
and measured H/3 equivalent widths (W(H/3)). We mea- 
sure W(H/3) to be ~100A for both regions, suggesting a 
burst age of ~4Myr dSchaerer fc Vaccalll99cl their Figure 



Table 6. O star populations and N(WR)/N(0) ratios for indi- 
vidual clusters from UV/optical spectroscopy using Starburst99 
spectral synthesis models, plus O star content from giant HII re- 
gions from He* imaging, after correction for ionizing fluxes from 
WR stars. We adopt N(LyC)=48.9 for LMC metallicity 07V 
stars, and estimate r\ from W(H/3) for regions A and B (see text) 



Region Al A2 


Bl + 2 


UV/optical Spectroscopy 


M(xlO 5 M ) 2.0 2.2 
N(O) 550 600 
N(WR)/N(Q) 0.2 0.1 


1.6 
450 
0.1 


Ha Imaging 


logQjf 6 52.39 
N(07V)t 2000 
N(0)t 4000 


52.19 
1600 
3200 



-0.5 



7). For a burst of this a ge and a Salpeter IMF, 77 is 
JSchaerer fc Vaccalll99cl their Figure 21). 

We present log Qo bs , N07V, No derived in this analysis 
in Table |S] These exceed the continuum derived cluster O 
star content by an order of magnitude. These would have 
been reduced by a factor of three had they been derived 
from Ha slit spectroscopy. These large difference suggest 
that NGC 3125 may host additional young massive clusters 
which are optically obscured. 

Recall from Figure 0^d) that region A is exceptionally 
bright at Ks. Our combined SB99 flux distributions for clus- 
ters Al and A2 predict that these clusters contribute only 
~30% of the observed Ks flux of region A. The addition of 
a visually obscured cluster (Eb-v ^ 1.5 mag), spatially co- 
incident with the optically visible clusters may resolve the 
observed IR excess for region A. To reconcile the high num- 
ber of equivalent 07V stars derived from the Ha image, 
a cluster with mass ~8x10 5 Mq and age ~l-2Myr would 
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Table 7. Comparison between o bserved He II 4686 and C IV 5808 
line flu xes in regions A and B bv lKunth fc Sargend Jl98lL K S811. 
KS81). IVacca fc Contil Il992t VC921. ISchaerer et al.l Il999t S99) 
and the present study, including differences in interstellar red- 
dening, with line luminosities calculated for a common distance 
of 11.5Mpc. 



KS81 VC92 S99 This study 

E fl _ v (Gal+Int) 0.4 0.08+0.4 0.25+0.27 0.08+0.16 

F(4686)a 8.2x10- 15 7.32xl0~ 15 8.2x10~ 15 8.5xl0~ 15 

L(4686)a 1.3x10 37 6.1x10 38 7.9x10 38 2.7x10 38 

F(5801)a 3.5x10" 15 2.0xl0~ 15 

L(5801)a 2.2x10 38 6.0x10 37 

E fl _v(Gal+Int) 0.08+0.64 0.25+0.44 0.08+0.13 

F(4686) fl 4.3 xlO -15 4.9xl0" 15 4.0xl0" 15 

L(4686)s 8.2xl0 38 8.4xl0 38 1.3xl0 38 

F(5801) s 6.0xl0~ 15 2.8xl0" 15 

L(5801) s 2.2xl0 38 7.6xl0 37 



be necessary. Indeed, visually obscured, young massive clus- 
ters are com mon in dwarf irregu lar starburst galaxies such 
as NG C 5253 jTurner et alJl200fJ) and He 2-10 IVacca et alJ 
120021) . 

For region B, clusters Bl and B2 are predicted to con- 
tribute ~60% of the observed Ks flux. We do not consider 
it likely that there is an additional obscured cluster since ei- 
ther Bl or B2 may be somewhat older than 4Myr, such that 
red supergiants would contribute to the IR excess. In this 
case, the IR bright sources to the NW of the UV/optically 
bright clusters would dominate the Ha ionization from re- 
gion B if each possess a mass of ~2xlO 5 M0 and an age of 
l-2Myr. 



5 COMPARISON WITH PREVIOUS STUDIES 

We will now compare our derived massive stellar content for 
NGC 3125 with those published in the lit erature. Optically 
derive d properties will b e com pared with Kunth & Sargent 
]l98ll) . IVacca fc Contil ]1992T) and ISchaerer et alJ il99Sf) 
whereas UV comparisons will be made to the HST/STIS 
survey of lChandar et alJ (1200411 . 

5.1 Optical Studies 

Table compares the measured line fluxes for He II A4686 
and CrvA5808 for regions A and B in the present study 
to those from the literature, showing, in general, very good 
agreement. Differences in line luminosities (for a uniform 
distance), and hence WR con tent, relate primar i ly to red- 
dening. Recall from Section |H| iKunth fc Sargentl lll98ll) de- 
rived reddenings from high Ba lmer lines, negl e cting cor- 
rections for stellar absorptions. IVacca fc Contil (I1992T) ob- 
tained red denings from Hq/H/3 observed during different 
conditions. ISchaerer et al] <ll999T) applied erroneous redden- 
ings from the literature, including a foreground extinction 
of E s _v=0.25 rather than Ay=3.1E s _y =0.25, although 
their assumed total extinction fortuitously closely agrees 
with IVacca fc Conti 

The major uncertainty in the derived WR content (for 



our derived reddening) result from the adopted line luminos - 
ity calibration. For WN subtypes. IKunth fc Sargentl ]l98lt) 
adopted a dominant mid-type WN population due to the 
marginal detection of Niv A4058 emission. We support this 
on the basis that both N IV AA4603-20 and N ill AA4634-41 
are weak due to the contribution b y C HI A4650 from WC 
star s (recall Figure [p. In contrast, IVacca fc Contil dl992l) 
and ISchaerer et al] il999T) adopted a dominant late-type 
WN population, albeit with a similar Hell line luminosity 
due to their inclusion of WN 6 subtypes. 

For WC subtypes, both ISchaerer et all lll999F) and the 
present study adopt a dominant early-type WC popula- 
tion on the basis that Cm A5696 is weak/absent. Both 
studies adopt a sim i lar C IV A5808 line luminosity as 
Crowther fc Hadfield 2006) support the earlier result of 
Smith et alJ ll990b1 using an increased LMC sample. 

Previous studies of the O star populations have 
been based on H / 3 flux es d erived from slit spect roscopy. 
IKunth fc Sargentl il98 Jl and lVacca fc Conti lll992ft assume 
WR stars do not contribute to the Lyman continuum and 
estima te O star populations of ~2000-2400. ISchaerer et al] 
( 1999) derived O star populations from H/3 spectroscopic 
line fluxes, but with the addition of correcting for the WR 
contribution and the evolution of the O star population, 
revising the O star numbers for NGC3125-A and -B to 
3240-6470 and 3450-6900, respectively. In contrast, this 
study has derived O star populations using Ha imaging 
for the Giant Hn regions, plus Starburst99 continuum fits 
to UV /optical spectroscopy for the clusters. For the latter, 
plus reduced reddenings, we estimate N(O) +2 = 1 1^0 an d 
N(0)B1 + 2 = 450- Differences between these and previously 
published results are again primarily attributed to adopting 
lower reddenings. 

5.2 UV Studies 

In Section 14.1.21 we have derived a content of 110 WN5- 
6 stars for NGC 3125-A1 from the HST/STIS Hell A1640 
line flux, whilst IChandar et al] derived a content of 5000 
late WN stars from the same dataset (using an extrac- 
tion window of 15 pixels versus 13 pixels here, and neglect- 
ing slit losses). The origin of this major difference is not 
due to the assumed intrinsic Hell A1640 line luminosity, 
since we assume a higher value o f 1.8 xlO 37 ergs -1 per WN 
st ar jCrowther fc _ Hadfieldl I2006T) versus 1.2xl0 37 adopted 
bv lCh^n^ar^t^lTil2^04T) ~As with the optical data, we have 
found that the difference bet ween WR popu lations derived 
here and those published by IChandar et al.l can be readily 
explained by differences in the adopted internal extinction 

law. 

The lCalzetti et al.l (Il994l) law was obtained using IUE 
large aperture (10x20") observations of distant starbursts 
with a median distance of 60 Mpc, i.e. a physical scale of 
~3 kpc. These observations sample a composite of stars and 
gas, sufferin g different extin c tion p roperties. In contrast, the 
SMC law of iBouchet et al.l dl985h was obtained from IUE 
observations of individual stars at a distance of ~60kpc, 
i.e. a physical scale of 3pc. Consequently, our HST/STIS 
spectroscopy (0.2x0.3") of cluster Al, sampling a physical 
scale of 10xl5pc, is most naturally suited to a low metal- 
licity stellar extinction law since the aperture is dominated 
by stellar extinction properties rather than interstellar gas. 
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Figure 8. Comparison between the complete slit-loss cor- 
rected HST/STIS spectrum of Al (faint solid) and 4Myr, LMC- 
mctal licity SB99 flu x distributions reddened according to (a) 
ICalzetti et al l Il994l) starbu r st ext inction law with Eg^ v =0.52 
(solidl: (b) iHowarth et alJ Jl983t) LMC extinction law with 
E™ T r =0.25 (dotted): (cl IBouchet et alJ Il985tl SMC extinction 
law with E™^^=0.16 (dashed). A Milky Way foreground extinc- 
tion of E(B-V)=0.08 has been applied in all cases and have been 
normalazied at 1600A. 



In Figure [S] we compare the complete STIS spectrum of 
cluster Al with various 4Myr, LMC metallicity SB99 SED 
(Section l4.2H that have been reddened to reproduce the ob- 
served 1250 - 1600A continuum flux distribution. For the 
first example the SED has been reddened using a starburst 
ext inction law for which we obtain E™_T V =0.52, as derived 
bv IChandar et al.1 (120041) . A cluster mass of 2xl0 6 M Q is 
needed to match the observed 1600A continuum flux. How- 
ever, the SED is highly discrepant in the near-UV, in the 
sense that the predict ed continuum fl ux is 60% too high. In- 
deed, if had selected a lChandar et alj had selected a different 
wavelength interval (e.g. 1600-2000A) they would have de- 
rived a signi ficantly lower reddening (e.g. E™2v=0.4mag) 
for the same ICalzetti et al] law. 

In Figure[S]we also include SB99 spectral energy distri- 
butio ns reddened to match the 1250-1 600A slope using the 
SMC (IBouchet et aljll98Eft and LMC iHowarth et al.lll983t) 
extinction laws. The former shows good overall agreement, 
using Eg™V=0.16 and M = 2.0 x 1O 5 M , including the weak 
2 175 A feature and gives the closest match to the observed 
continuum flux in the near-UV. For the LMC law, the far- 
UV slope is reproduced with E^^O.25 (3.5 xl0 5 M Q ). In 
this case, however, the 2175A feature is overestimated and 
the near-UV continuum flux is poorer than that for the SMC 
extinction law. 

Consequently, it is not possible to derive a consistent 
reddening for NGC3125-A1 using a starburst SED with ei- 
ther a LMC or Calzetti extinction law. In contrast, we obtain 
an excellent match to the observed UV SED for our optically 
derived reddening with a starburst spectral energy distribu- 
tion and a SMC extinction law and an internal extinction in 
excellent agreement with our nebular derived Ha/H/3 red- 
dening. 



6 CONCLUSIONS 

We have demonstrated that the WR populations of the two 
regions NGC 3125-A and B are substantially lower than pre- 
vious optical studies, and dramatically lower than previ- 
ous UV studies of the bright cluster Al within the giant 
Hn region A. Indeed, previous highly discrepant UV and 
optical results for NGC 3125-A may be reconciled using a 
Ha:H/3 derived reddening and a SMC extinction law for the 
internal extinction. We have obtained refined WR popula- 
tions in these regions by applying template spectra of typical 
WN and WC stars to the blue A4686 bump, which indicate 
that WC stars, may contribute significantly to the observed 
A4640 feature, commonly attributed solely to N ill A4640 
from late WN stars. 

IChandar et al] i2004T) argued that the presence of 
strong Hell A1640 emission in Al implies an exceptional 
WR population with N(WR)/N(0)>1. However, evolution- 
ary models for single stars at LMC metallicities predict a 
maxi mum N(WR)/N(Q) ratio of ~0.1 for an instantaneous 
burst (Sc haerer fc Vaccal fl998). whilst we estimate values of 
N(WR)/N(O)~0. 1-0.2 for clusters within the giant Hn re- 
gions NGC 3125-A and B. Consequently, our results broadly 
reconcile the observed massive stellar content of NGC 3125- 
A and B with evolutionary predictions for a young LMC 
metallicity starburst. 

The primary difference between the two UV studies was 
the choice of extinction law. A standard starburst extinction 
law is ideally suited to spatially unresolved stellar galaxies, 
such as high z LBGs, but when combined with a starburst 
spectral energy distribution is unable to reproduce the UV 
spectru m of a resolved star cluster such as NGC 3125-A1 or 
Tol89-l JSidoli et al]l2006h . In such cases, a SMC (or LMC) 
reddening law is necessary to correct the complete UV spec- 
trum for extinction. Consequently, Hell A1640 resu lts for 
other nearby WR clusters studied bv lChandar et all should 
be treated with caution for the present. 

We have also demonstrated that the O star content de- 
rived from Hcv narrow band imaging is substantially higher 
than that estimated from continuum flux techniques. On the 
basis of near-IR imaging of NGC 3125 we propose that the 
discrepancy for region A can be resolved by the presence of 
an additional cluster which is optically obscured. Such clus- 
ters appear to be common in LMC- like metallicity starburst 
galaxies such as N GC 5253 and He2-10 jTurner et al]l200d : 
IVacca et al]l2002ft . For region B, the O star content may be 
resolved by the identification of two additional knots to the 
NW of region B which are bright at near-IR wavelengths. 

The present study emphasises the need for the highest 
possible spatial resolution. We have attempted to evaluate 
the properties of the two optically visible clusters within re- 
gion A based upon a pair of broadband ultraviolet (FOC) 
and far-red (STIS) flux ratios, plus the properties of the 
clusters within region B from a single, composite ultra- 
violet measurement. Refined properties for clusters within 
NGC 3125 require high spatial resolution datasets spanning 
the ultraviolet to near-infrared, as recently undertaken by 
program GO 10400 (R. Chandar, P.I.) using the HST Ad- 
vanced Camera for Surveys (ACS) and Near infrared Cam- 
era and Multi-Object Spectrometer (NICMOS). 
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